Fluctuations and redundancy in optimal transport networks 
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The structure of networks that provide optimal transport properties has been investigated in a 
variety of contexts. While many different formulations of this problem have been considered, it is 
recurrently found that optimal networks are trees. It is shown here that this result is contingent on 
the assumption of a stationary flow through the network. When time variations or fluctuations are 
allowed for, a different class of optimal structures is found, which share the hierarchical organization 
of trees yet contain loops. The transitions between different network topologies as the parameters 
of the problem vary are examined. These results may have strong implications for the structure and 
formation of natural networks, as is illustrated by the example of leaf venation networks. 

PACS numbers: 89.75.Hc, 89.75.Fb, 05.65.+b, 89.75.Da 
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The structure of networks that provide optimal trans- 
port efficiency is a question of long-standing interest [ll, 
S S| , with applications ranging from the design of infras- 
tructures such as water and power supply networks to the 
analysis of natural networks such as the vascular system 
of plants [1] and animals ^ or river basins Q. Optimal 
networks may be analyzed in terms of their topology, 
and, if they are large enough, in terms of their statistical 
properties. In this respect, the self-similar structure ob- 
served in many natural networks has strongly stimulated 
the general investigation of optimal networks. For nat- 
ural networks, the discussion of their functional proper- 
ties is also intimately related with that of the patterning 
mechanisms by which they arise. 

While different formulations of the optimal transport 
problem have been developed in different contexts, they 
recurrently yield tree-like structures, with a single path 
connecting any two points 

BBSS 0- Many nat- 
ural networks conform with this predicted tree-like or- 
ganization. However, it stands in strong contrast to the 
structure of leaf venation networks, which contain many 
loops jTlJ . This structure could be interpreted as a com- 
promise between transport efficiency and other require- 
ments, such as tolerance to damage. On the other hand, 
it has been suggested that the redundancy of leaf ve- 
nation networks could be an adaptation to the varying 
water demands of different parts of the leaf [l3| . In this 
Letter, we examine the structure of optimal networks car- 
rying a flow that varies with time or fluctuates. Beyond 
the particular example that motivated this study, this is 
a question of general relevance, since the stationary flow 
usually assumed in the study of optimal networks is al- 
ways an idealization. We show that fluctuations can give 
rise to a different class of optimal structures, which con- 
tain loops while retaining the hierarchical organization of 
trees. Introducing a suitable measure of network redun- 
dancy, we characterize the transitions between different 
network topologies as the parameters of the problem are 
varied. We discuss the implications of these flndings for 



the structure and formation of natural networks, showing 
how they can account for the presence of loops in some 
cases and their absence in others. 

When considering networks that transport fluid or 
electrical current, a natural formulation of the optimiza- 
tion problem consists in minimizing the total dissipation 
rate with a limited amount of resources [1, Q . Consider 
a network formed of vertices k and edges (fc, I) having 
conductances Hki- For simplicity, we consider a regular 
network (all edges have length one). Sources ik are con- 
nected to the vertices, with J^k ~ and the edges 
carry currents Iki , which by Kirchhoff 's current law sat- 
isfy ik = X); ^ki- The problem consists in minimizing the 
dissipation 
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under the constraint that 
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The constant K'^ can be interpreted as the amount of ma- 
terial available to build the network, with the exponent 
7 depending on the nature of the network. For instance, 
7 = 1 for electrical wires or porous pipes. In his analy- 
sis of vascular networks, Murray Q considered Poiseuille 
flow in hollow pipes and a fixed volume of fluid, which 
corresponds to 7 = 1/2. 7 < 1 means that for a given 
cost, the dissipation can be lowered by concentrating the 
flow along links of large conductance. The case 7 > 1 
may be considered of little practical relevance, as it is 
then more economical to build several parallel links hav- 
ing a small conductance than a larger one of equivalent 
capacity. When the flow through the network is station- 
ary, it is found that for 7 < 1, all local minima of the 
dissipation are spanning trees, while for 7 > 1, there is 
a single, global minimum, and all edges have non-zero 
conductance 0, S] ■ 



2 



To study the influence of fluctuations, we now consider 
the state of the sources i = (ii, . . . ,iN) as a random 
variable with probabihty density p{i). The dissipation in 
state i is 



J(i) 
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where the Iki{i) are the currents in that state. The prob- 
lem is now to minimize the average dissipation 



< J >= 



J(i)dp(i). 



(4) 



Introducing a Lagrange multiplier A, the solution is ob- 
tained by minimizing 



^ =< J > -xy^ K 
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In many cases, the currents derive from a potential, such 
as pressure or electrical potential, and are uniquely deter- 
mined by the sources and conductances. However, this 
requirement need not be explicitly stated in the prob- 
lem. Indeed, following Thomson's principle 13], when 
arbitrary currents satisfying Kirchhoff's current law are 
considered, the currents that minimize the dissipation 
derive from a potential. In all cases, the currents and 
conductances may thus be varied independently in the 
optimization. Minimizing S with respect to k^i yields 



< Hi > 
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This is our first important result, which relates the con- 
ductance of each link in an optimal network to the mean 
square current that it carries, generalizing a similar rela- 
tion derived for stationary flows [8|. 

The above relation also serves as the basis for our 
numerical investigation of the structure of optimal net- 
works. Indeed, local minima of the dissipation can be 
computed by iterating Eq. [5] _8|. For a given set of 
conductances, the < /^^ > can be determined from the 
source correlations. The potentials Ufc(i) are related to i 

by 



or in matrix notation i = Cu(i), and 

< UkUi >~ RkmRln < irnin >, 

where R = C^^. Finally, 

< 4 >= 4i {<ul> + <uf > -2< UkUi >) 

The above minimization procedure is started with ran- 
dom conductances. 
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In what follows, an n x n square network {N = n^), 
with — 1 sources (2 < fc < N) and a sink [k = 1) in 
one corner is considered. The sources have unit average 
and uncorrelated fluctuations of amplitude a: < ik >= 1 
and < ikii >= 1 + cr^Ski, 2 < A:,Z < A^, where S de- 
notes the Kronecker delta function (the < ikii > are 
obtained from Efc*fe — 0)- For sources that are ran- 
domly switched on and off with equal probability, a — 1. 
The network topologies observed with constant sources, 
trees and rather uniform networks with many loops, are 
recovered for 7 small enough (Fig. [Ha)) and for 7 > 1 
(Fig. [T]^c)), respectively. Our central result is that in an 
intermediate range of 7, however, hierarchical networks 
with loops are obtained (Fig. [Hb)). An intuitive inter- 
pretation of this result is that 7 < 1 favors grouping the 
flow along edges of high conductance, while loops provide 
alternate routes to accommodate fluctuations. The local 
minima obtained for 7 < 1 depend on the initial condi- 
tions. For small systems, however, we have checked that 
local minima containing loops yield a lower dissipation 
than the optimal tree by constructing all spanning trees 
rooted at the sink [l^. Note also that local minima are 
relevant where self-organized processes are concerned, as 
in the case of natural networks. 

To investigate the transitions between these differ- 
ent network topologies, it is suitable to introduce a 
measure of network redundancy. A simple such mea- 
sure is the number of loops L, which can be deter- 
mined using Euler's formula, L — I + number of edges — 
number of vertices, where only edges that exceed a 
threshold conductance are counted. Fig. [2t^a) shows 
that loops appear above a finite threshold 7c, and their 
number progressively increases to reach its maximum as 
7^1. The number of loops, however, is not an ideal 
measure of network redundancy, as it is a purely topolog- 
ical one, and depends on an arbitrary threshold. These 
limitations can be overcome by introducing an entropy- 
based measure of redundancy, which is similar to mea- 
sures of reliability appearing in the literature on water 
distribution networks [l5l |. and could be of general in- 
terest in the analysis of complex networks. First, the 
redundancy rat between two vertices a and b is defined, 
using an analogy with random walks. Consider the cur- 
rents Iki corresponding to a unit current injected into 
the network at a and extracted at b {ik — 5ka ~ 5kb) and 
transition probabilities 



Pkl^l ^'.^>o'''"' 

[ 0, hi < 0. 



(10) 



The probability of any path k = (fci, . . . ,ki) from a to 6 
{ki — a, ki = b) is 



p(k)= n Pk^k,^,, 
i<j<i 



(11) 



and the redundancy Tab is defined as the entropy of the 
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FIG. 1: Network structures obtained with fluctuating sources for different values of the parameter 7 (system size n = 16, 
fluctuation amplitude cr = 1). In each network, the sink is at the lower left corner. The width of each edge is proportional to 
the square root of its conductance (corresponding to the diameter of electrical wires or porous pipes), (a) Tree-like network 
(7 — .25). (b) Hierarchical network with loops (7 = .75). (c) Network with many loops and no hierarchical organization 
(7 = 1.25; notice that the network is not perfectly uniform, conductances being higher near the sink). 



path distribution: 



(12) 



This expression can be transformed into a sum over the 
nodes, r^fc = J^i^i^i^ where Pi is the total inflow into 
node i and Si = — J^j Pij Pij ■ '''ab is equal to if 
there is a single path from a to 6, to log n if there are n 
independent paths of equal conductance from a to 6, and 
can take arbitrary positive values in general situations 
(e''"'' can be interpreted as a generalized number of paths 
between a and h) . In the present case of identical sources 
and a sink, a sensible measure of the overall redundancy 
R of the network is the average redundancy between any 
source and the sink, i.e. 



R 
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As shown on figure [D^b) , the variation of the redun- 
dancy R with the parameter 7 is qualitatively similar 
to that of the number or loops. The redundancy is 
when there are no loops, and gradually increases above 
the threshold 7c. Note however that it is not constant 
for 7 > 1, and only progressively tends to its maximum 
as 7 — > +00 and the conductances become uniform. The 
value of the maximum is of order n, reflecting the ex- 
ponential growth of the number of possible paths in the 
network. 

Using this measure of redundancy, we now turn to 
an analysis of the transitions between different network 
structures as the parameters of the model are varied. 
Fig. [3] shows the dependence of the redundancy R on the 
amplitude of the fluctuations. Expectedly, the transition 
is sharper (7c is larger) for smaller amplitudes. Although 
a limited range is accessible numerically (each step of 
the minimization procedure involves the inversion of an 
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FIG. 2: Evolution of network structure vs. 7 with constant 
(crosses) and fluctuating (lines) sources {n — 8, a — 1, av- 
erages over 100 realizations), (a) Number of loops L. For 
7 > 1, all conductances are non-zero and L = {n-lf. (b) 
Redundancy R. Note that R is not constant for 7 > 1. 
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FIG. 3: Redundancy R vs. 7 for different fluctuation ampli- 
tudes (n = 8). 



Nx N matrix), we have also examined the dependence of 
optimal network structure on system size. Fig. [4] shows 
the number of loops and redundancy as a function of 7 
for different system sizes. We find that 7c appears to re- 
main constant. The number of loops scales as (n — 1)^, 
suggesting a constant average loop area. On the other 
hand, the normalized redundancy appears to increase less 
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FIG. 4: Dependence of network structure on system size 
(a — 1). (a) Number of loops L rescaled by l/(n — 1)^. (b) 
Redundancy R rescaled by 1/n. 



rapidly above 7c for larger values of n. This may be due 
to a weaker redundancy on the large scales, due to the 
averaging out of the uncorrelated fluctuations. 

In summary, we have shown that fluctuations dramati- 
cally change the optimal structure of transport networks, 
favoring the presence of loops. The optimal structure is 
determined by the correlations between the currents in- 
jected into the network. We restricted our analysis to un- 
correlated fluctuations, and further investigations should 
address the influence of correlations, which can be ex- 
pected to increase redundancy on larger scales. It would 
also be of interest to examine the statistical properties 
of the structures obtained. One the one hand, optimal 
trees are generally found to be self-similar Q . This prop- 
erty might allow an estimation of the threshold for the 
appearance of loops by balancing the cost of inserting a 
new link and the corresponding dissipation reduction. On 
the other hand, one could ask whether self-similar struc- 
tures can also be obtained in reticulate networks. This 
might be the case when the fluctuations exhibit power- 
law correlations. In another line, it has been shown for 
stationary flows that the geometry of network junctions 
satisfies special relations when dissipation is minimized 
with respect to the locations of the vertices fl6l| , and one 
could attempt to generalize these relations in the pres- 
ence of fluctuations. 

An application of our results may be found in the struc- 
ture of natural networks, and in particular in the retic- 
ulate structure of leaf venation networks. Indeed, the 
rate of transpiration in the different regions of a leaf 
is known to fluctuate over time 17|. In addition, leaf 



veins are bundles of smaller vessels, with transport prop- 
erties similar to a porous medium, corresponding to a 
value of 7 close to 1 (yet smaller than 1 since the ves- 
sel diameters increase with vein diameter). This is the 
range where the structure of optimal networks is most 
sensitive to fluctuations. Interestingly, long-range cor- 
relations have been observed in the fluctuations of leaf 
transpiration rate ^is'l, which could explain why leaf ve- 
nation networks exhibit redundancy at every scale. On 
the other hand, the vascular networks of animals can be 
described by the smaller value 7 — 1/2, and a tree-like 



structure would be predicted even if moderate fluctua- 
tions are taken into account. Our results may also be 
relevant to the formation of reticulate networks such as 
leaf venation networks. Indeed, Eq. [6]shows that optimal 
structures may be obtained through local adaptation of 
the conductance to the flow. Returning to a functional 
perspective, another benefit of redundancy is tolerance 
to damage, and it would be of interest to compare the 
structures obtained here with those found when explic- 
itly optimizing for tolerance to damage. 
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